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Abstract 

A self-consistent description of the MIT bag parameters as functions of the nu- 
clear matter density is presented. The subnuclear degrees of freedom are treated in 
the Quark-Meson Coupling Model, considering the equilibrium conditions for the 
bag in the nuclear medium. The hadronic interaction is described in the framework 
of the quantum field theory of hadrons through several models. We have obtained 
the behavior of the bag radius and the bag parameters B and Zo, taking their deriva- 
tives with respect to the mean field value of the scalar meson as free parameters. A 
discussion on the variation range of these derivatives is given. 
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It is well known that Quantum Chromodynamics (QCD) is the most accepted candi- 
date for a theory of strong interactions. Although QCD has been successful in describing 
the high energy domain, its application seems difficult in the energy range associated with 
nuclear phenomena. Attempts to sort this difficulty have been carried out by using effec- 
tive models, among them the so-called bag models are widely used [J. They provide 
an acceptable description of the free hadronic properties. The energy-momentum conser- 
vation of an isolated bag is imposed through the so-called non-linear boundary condition 

Si- 

The Quark-Meson Coupling Model (QMC) has been developed since the early work 
of Guichon ||. In this model the nucleons are described as non-overlapping MIT bags 
confining the quarks inside them. The quarks interact through the exchange of scalar (a) 
and vector (u^) mesons. Several refinements and applications of this model have been 
done for instance the QMC model has been used to describe nuclear and neutronic 

matter j|, |5|, as well as finite nuclei || ||. 

On the other hand, the nuclear phenomena has been treated in the context of quan- 
tum field theory by using structureless hadronic fields as the effective degrees of freedom. 
This subject reached interest since the pionnering work of Walecka ||. Different aspects 
of the phenomenology of nuclear matter and finite nuclei have been studied in the frame- 
work of the so-called Quantum Hadrodynamics (QHD) with successful results 0. The 
original Walecka model has been extended with additional polynomic potentials |jT0[ and 
alternative non-polynomic interactions in the scalar channel |ll||-[fL4||. 

If the validity of the equation of state , effective mass, etc., derived from an effective 
hadronic lagrangian is assumed, one can ask about the hadronic substructure consistent 
with those properties. We are interested in the implications that the QHD description have 
for a picture which deals with subnuclear degrees of freedom. A clear relationship between 
the QMC model and the Walecka model has been stablished by Saito and Thomas (3| . 
In the present work we have studied the reliability of a coherent description of in-medium 
nucleon properties in terms of QHD and QMC. 

We briefly recall some basic features of the QMC model |3|, |J. If isospin symme- 
try breaking is neglected the meson fields cr(x) and cu M (x) are sufficient to describe the 
problem. Inside the bag the equation of motion for quarks of mass m q is given by 

(i0 - m q )m q (x) = l-g q Mx) + g^{x)]^ q {x), (I ) 

where g% and are the quark-meson coupling constants associated with the a and 
fields, respectively. In the Mean Field Approximation (MFA) the meson fields are replaced 
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by their mean values, which become constants in infinite nuclear matter, i.e. a = a and 

The normalized quark wave function for a spherical bag of radius R is 

».(f, t )=^-xf ^ ) )**=, (2) 

y iPga-rjiiyr/R) J V4vr 

where r = \f\, Xq is the quark spinor and the normalization constant M is given in ||. 
We have introduced the effective quark mass m* = m q — g^a and the effective energy 



eigenvalue e q = VL/R + g^co, where f2 = yy 2 + (Rm*) 2 . The y variable is fixed by the 
boundary condition at the bag surface jo{y) = (3 q ji(y) as in reference [TJ and f3 q = 
y/(n-Rm*)/(n + Rm*). 
The bag energy is given by 

E b = + \itBR* (3) 

where B is the energy per unit of volume and Zq takes into account the zero point energy 
of the bag. The nucleon mass is defined by including the correction due to the spurious 
center of mass motion [fTfl 



Ml = ^El - 3(y/RY. (4) 

It is usual to determine the bag parameters at zero baryon density to reproduce the 
experimental nucleon mass Mb = 939 MeV. Simultaneously it is required the equilibrium 
condition for the bag dM}>(a)/dR = 0. In [EJ the bag parameters B and Zq are constants, 
although it was found a relative change of 1% in the radius at the saturation density 
(p = 0.15 fm -3 ) as compared with its value at zero baryon density. 

In the normal QMC model the bags interact by the same mechanism that couples 
mesons to quarks inside a bag. In the present work we have proposed a description of 
the hadronic interaction by a set of different effective models commonly used in QHD. 
In all the cases considered we have used a linear coupling between the nucleon field 
ty(x) and the vector-meson field, with strength g w . However the models differ in the 
interaction term ^(x)VN a ^(x), between the nucleon and the scalar meson field. In Table 
I we show the explicit form of for these four models. The coupling constants g a and 
g^ have been fixed to reproduce the binding energy per nucleon, E b = 16 MeV, and the 
nuclear saturation density, p = 0.15 fm -3 , in the MFA. The isothermal compressibility 
k = 9p (dPh/dpB)r is a usefull quantity in order to analyze the fitness of hadronic models, 
since its value at the nuclear saturation density has been well determined to range between 
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100 — 300 MeV [15 . The isothermal compressibility evaluated in the MFA is also shown 
in Table I ill. 



The Euler-Lagrange equations for nucleons and mesons in the MFA are given by 

{i$ - M N )*(x) = (gjsyo - V Na )^{x), (5) 

m a a = —^-{e = °) Ps, (6) 
mlo = g^p, (7) 

where p s =< 1 l r (x)\E'(x) >, p is the baryon density; M N , m a and are the masses of the 
free nucleon, the scalar and vector mesons, respectively. 
At zero temperature we have 

' - §• (8 > 

and 

Ps = 7K-rJ d 3 kQ(k F - \k\) » , (9) 
( 2?r ) J \ M* N 2 + fc 2 



where kp is the nucleon Fermi momentum. The effective nucleon mass and the energy 
spectrum e(k) are given by 

M* N = M N - V Na , (10) 

and 



e{k) = ^/Mf + k^ + g^. (11) 

We have used different notations for the nucleon mass entering in QHD (Mjy-), and the 
nucleon mass generated by the bag model (M£). 

Equation (g) is a self-consistent definition for the mean field value a, indeed from this 
equation we see that the derivative dM* N jda determines the dynamics of the scalar field. 
The QMC and QHD descriptions produce coherent results if the following equation is 
fulfilled 

M* N {a) = M;(a), (12) 

together with g a = 3g q a , g u = 3g«, |fl. 

The stability of the bag in the nuclear medium with respect to volume changes is 
imposed by 

P b {a) = P h (a), (13) 

where Pb{&) is the internal pressure generated by the quark dynamics and Ph{°') is t ne 
external hadronic pressure. 
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The bag pressure can be written as 



E b S t , 16 d2d 3m^\ 3y 



+ ^R'B + -^-\- -r^TTZ, (14) 



4ttR 2 M* \ R 3 tt J 4ttR 5 M* 



and the pressure for uniform nuclear matter is given by P h = -\T % \ where T" is the 
trace over the spatial components of the energy-momentum tensor 



2M* N 3 

-Ol — 7T- 

7T 



kp (5e F e| \ , ( e F + kj 



-In 



M* N \ 24 YlM^j 8 "" V M* N 



(15) 



The equation ( O ) is a statistical equilibrium condition on the bag surface which 
ensures a direct relation between nuclear matter bulk properties and the stability of the 
confining volume. 



The equations ([12]) and (13) can be used to obtain the bag parameters as functions 
of the density. However, they are not enough to determine all the functions involved. 
Henceforth we have included the additional equations obtained from the derivatives of 
the boundary condition and of the equations (|l|) and (|13|), taking y,m*, R, z and B 
as independent functions. In fact we have considered A = dBjda and fi = dzo/da 
as constants. This procedure is equivalent to take only the linear contributions in the 
expansion of the functional relations ( |i~2"D and (|13|). 

After the hadronic coupling constants have been adjusted to reproduce the saturation 
conditions for nuclear matter, equation (^) can be used to obtain the solution a for each 
baryon density. Equations @, (f|), ( |i~0D and (0) can be used together with eqs. flllf ), ([14]) 
and (|T5| ) to determine B and Zq when the value of the in-medium bag radius is provided. 
The last quantity is obtained by solving self-consistently the additional equations for fixed 
values of the parameters A and /i. To search for appropriate values of A and /i we have 
evaluated R, B and z at zero baryon density as a function of (A,/i) using the model 1. 
A drastic change in the behavior of the quantities considered is found when fi goes from 
negative to positive values. 

The bag radius R at the saturation density has been studied as a function of A. For 
the hadronic models considered here and taking the bag radius at zero baryon density 
Rq = 0.6 fm we have found that R/Ro > 1 only for a restricted range of fi. For the 
following discussion we have taken two sets of values; set I (A = —5.28 fm~ 3 , \i = —0.50 
fm) and set II (A = 0,/i = 1.6 fm). 

The bag radius as a function of the density is shown in fig. 1, the quark mass at 
zero baryon density has been fixed at m q = 10 MeV. It can be seen that set I gives an 
asymptotical constant bag radius for every model, whose values do not depend on the 
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details of the interaction and it is diminished as compared with its vacuum value. The 
set II provides a model dependent radius at high densities, in this case models 1 and 2 
predict a breaking-down of the bag picture. Models 3 and 4 have a stable behavior for all 
the densities considered here, even when set II is used. 

In figs. 2 and 3 we present the density dependence of _B 1//4 and of z , respectively. 
Models 1 and 2 exhibit an opposite behavior for B 1 ^, when set I or set II are used. 

For model 1 and set I B takes negative values as the baryon density is sufficiently in- 
creased, thus the bag bulk energy must decrease with increasing volume. The increment 
of the kinetic energy compensates this fact, giving a slowly decreasing total bag energy 
and a stable bag radius (see fig. 1). On the other hand, B grows drastically at high 
densities for set II, a small volume increment gives rise to a large increment in the bulk 
energy. Therefore in order to get a slowly decreasing M b *, the bag radius R must decrease 
at the same rate as 

£1/3 

grows. When R approaches to zero a subtle cancellation among 
the quark kinetic energy, the zero point motion parameter and the center of mass cor- 
rection takes place. The raising of Zq (see fig. 3) is not sufficient to reach the dynamical 
equilibrium and hence the system reduces its volume as far as possible. 

The steep behavior of models 1 and 2 as compared with models 3 and 4 is due to the 
fact that the first mentioned models give a stiff equation of state and a fast decrease for 
the effective nucleon mass. 



Our results can be compared with those obtained by Jin and Jennings In their 
work a phenomenological description of B as a function of the density is given in terms of 
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a set of two free parameters with no direct dynamical interpretation, namely 5 and B, 
for model I of [0 . Starting from the usual approach in the QMC model, in ref . [[l|J the 
constant B is replaced by a function of the baryon density. Within this framework they 
found a monotonous density dependence of the bag parameter B. For 5 = 4 and = 1, 
the ratio B/B decreases 80% at the nuclear saturation density and the corresponding 
in-medium bag radius increases 60%. In our approach we have obtained an increase of 
12% for the bag radius (fig. 1) and a corresponding decrease of 19% in B (fig. 2), when 
model 4 and the set II of parameters are used. 

In this work we have studied the coherence of the QHD and QMC descriptions by 
using the equilibrium conditions for the bag in nuclear matter. Thus we have stablished a 
link between these frameworks which is able to explain the nucleon substructure changes 
taking into account the nuclear matter equation of state, the nucleon effective mass, etc. 

We have evaluated the density dependence of the bag parameters and the bag radius. 
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In our model we use two dynamical quantities, i.e. the derivatives dB/da and dz^/da as 
free parameters and we have explored their possible variation range. We have found two 
different dynamical regimes for these parameters. 

The inclusion of thermal effects and the study of the EMC effect in this framework 



will be reported in future works fl? | 
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Table I. 



Model 


v Na 


9o 


9u 


K 


[MeV] 


1 




11.04 


13.74 


554 


2 


MNg a tanh(g a a / Mn) 


9.15 


10.52 


410 


3 


M N [l-exp(g a a/M N )} 


8.34 


8.19 


267 


4 


g a a/(l + g a a/M N ) 


7.84 


6.67 


224 



Table I: Nucleon-scalar meson interaction terms, coupling constants and isothermal 
compressibility k, for several effective hadronic models used in this work. 
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Figure Captions 



Figure 1: The bag radius R as a function of the relative baryon density p/po, where po 
is the nuclear matter density at saturation, for hadronic models indicated as M 1 
(Walecka model), M 2, M 3, and M 4 (Zimanyi-Moszkowski model). Full and dashed 
lines correspond to set I and set II of parameters, respectively. The quark mass at 
zero baryon density has been taken as m q = 10 MeV. 

Figure 2: B 1 ^ in [MeV] function of the relative baryon density p/po- The conven- 
tions and parameters used are the same as in Fig.l. 

Figure 3: z as a function of the relative baryon density p/po- The conventions and 
parameters used are the same as in Fig.l. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9611028v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9611028v2 



